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Abstract 

We consider the processes d + d n + ^He, d + ^He — > p + ^He, 
d + ^H ^ n + ^He, ^He + ^He ^ p + p + '^He, + '-^He ^ d + ^He, with 
particular attention for applications in fusion reactors. After a model in- 
dependent parametrization of the spin structure of the matrix elements for 
these processes at thermal colliding energies, in terms of partial amplitudes, 
we study polarization phenomena in the framework of a formalism of helicity 
amplitudes. The strong angular dependence of the final nuclei and of the 
polarization observables on the polarizations of the fuel components can be 
helpful in the design of the reactor shielding, blanket arrangement etc.. We 
analyze also the angular dependence of the neutron polarization for the pro- 
cesses d + d ^ n + ^He and d + ^H ^ n + "^He. 
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I. INTRODUCTION 



Nuclear fusion reactions, like d+d n + ^He, or d+^H — > n + ^He, are characterized by a 
large dependence on the spins of the colliding particles. It has been suggested |]T| to use this 
property in magnetic fusion reactors with polarized nuclear fuel. A magnetic field of about 
1 kG can keep the necessary direction of the polarization of the interacting nuclei, during 
a time which is longer in comparison with the reaction time. Different technical solutions 
might be used: injection of polarized frozen pellets, or polarized targets for inertial fusion. 

The strong dependence of the fusion reaction rates on the polarization states results in 
an increasing or a decreasing of the cross section (with respect to the unpolarized case), 
depending on the colliding nuclei polarization directions. These characteristics can be used 
to optimize a fusion reactor in different ways: 

• The possible enhancement of the fusion rates for d + ^He and the suppression of d + d- 
collisions would make this fuel competitive with d + ^if, as it would, in particular, 
result in a dean reactor. 

• The strong anistropy of the neutron angular dependence in ci + ^if-coUision helps in 
optimizing the reactor shielding and the blanket design. 

• d + ^iJ-collisions can be source of intensive monochromatic polarized neutrons, with 
the choice of the polarization direction. 

A precise knowledge of the spin structure of the threshold matrix elements for the pro- 
cesses induced by: d + ^He, d + ^H, ^He + ^He, +^ife-collisions is required. At energies 
up to 10 keV, which are typical for fusion reactors, the S-state interaction of the collid- 
ing particles dominates and the general analysis of polarization phenomena is essentially 
simplified. 

For this particular low energy domain, the general analysis of polarization phenomena 
(in terms of scalar or helicity amplitudes) does not seem adequate, as it does not include 
the simplified characteristics of threshold regime. The main reason is that, at threshold, 
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there is only one physical direction: the 3-momentum of the produced particles. Therefore, 
the parametrization of the corresponding matrix elements has to be done on the basis of 
one direction. It can not be derived as a limiting case of a general approach based on two 
independent directions, the 3-momenta of the initial and final particles, which define the 
reaction plane. For S-state collisions we can not define this plane, having only one direction. 
The analysis of threshold polarization phenomena is also simpler and requires a dedicated 
parametrization, which is not a limiting case of a general parametrization. 

Note that, in this connection, the spin structure for the threshold kinematics is equivalent 
to (and sometimes simpler than) the case of collinear kinematics. 

Our aim is to analyze here in the most general and complete form the reactions relevant 
to magnetic fusion reactors, with polarized fuel. Following the line of our previous paper 
0, where we considered the reaction d + ^He p + '^He, we will give here the general 
parametrization of the threshold amplitudes for d + ^He, d + d and ^He + ^i7e-collisions, with 
special attention to the angular distribution of the reaction products for different possible 
polarization states of the colliding particles, without any particular assumption about the 
reaction mechanism. For this aim we develop a formalism for the parametrization of the spin 
structure of the threshold matrix elements in terms of the small number of allowed partial 
amplitudes, taking into account the most general properties of the strong interaction. 

In a fusion reactor the reaction rates and the angular distributions depend on the direc- 
tion of the magnetic field. We use in this analysis a particular set of helicity amplitudes, 
with quantization axis along the direction of the magnetic field. (The effects of the magnetic 
field were not discussed in our previous paper 0). 

The paper is organized as follows. In section 2 on the basis of some of our previous 
results concerning the reaction d + ^H{^He) — > n{p) + '^He, we derive the angular dependence 
of the differential cross sections for different polarization states of the colliding particles, 
and the angular dependence of the polarization of the produced neutrons (protons). 

The process d + d ^ n + ^He, discussed in Section 3, is characterized by a set of three 
independent partial amplitudes (for S-state (irf-collisions) . We derive the limits of the integral 
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coefficients for polarized particles collisions, which quantify the change of the cross section 
due to the polarizations of the colliding particles. These coefficients depend only on the 
ratio of the square of the partial amplitudes. 

Section 4 contains the discussion of the properties of some processes induced by ^He+^He, 
+ ^H, and + '^ife-collisions. The Pauli principle (for colliding or produced identical 
fermions), essentially simplifies the spin structure of these reactions, and it is possible in 
some cases to give model-independent predictions for polarization phenomena. 

II. THE COMPLETE EXPERIMENT FOR THE REACTION 

A. Introductory remarks 

The reaction d n +^He in the near threshold region is very interesting for the 

production of thermonuclear energy and plays an important role in primordial nucleosynthe- 

3 + 

sis. The lowest - level of ^He has excitation energy = 16.75 MeV (only 50 keV above 
d +^if-threshold) and has a width of 76 keV. 

The microscopic explanation of the nature and the properties of this resonance is very 
complicated and still under debate in the physics of light nuclei. The interpretation |^ of 
this resonance as a shadow pole introduces a new concept in nuclear physics, after atomic 
and particle physics. The possibility that the corresponding shadow poles for the two charge 
symmetric systems d +^He and d + ^H (or p + '^He and n + '^He) occupy different Riemann 
sheets, due to the difference in electric charges of the participating particles, can not be 
presently ruled out. Such phenomena can be considered as a new mechanism of violation of 
isotopic invariance of the strong interaction [§. 

Due to the close connection of the three processes d + ^He ^ d + ^He, n + '^He ^ n + '^He 
and d + ^H ^ n + '^He through the unitarity condition, the partial wave analysis can not 
be performed independently for each reaction. The corresponding amplitudes are complex 



4 



functions of the excitation energy. The multilevel 7^— matrix approach allows to parametrize 

this dependence in terms of few parameters as shift, penetration factors and hard-sphere 

3 + 

phase shift 0. All characteristics of the = - -resonance, like the position, the width 



and particularly the Rieman sheet, can be found using an S*— matrix approach [^,|I^-[T2] 



The polarization phenomena are very important in the near threshold region, even for 
the S-state interaction. In this respect the reaction d + ^ n + ^He plays a special role, 
because the presence of a D-wave in the final state results in nonzero one-spin polarization 
observables, such as, for example, the tensor analyzing power. In order to fully determine 
the two possible threshold (complex) amplitudes, two-spin polarization observables have to 
be measured, for example in collisions of polarized deuteron with polarized ^i/e-target. Here 
we will generalize our previous analysis 0, taking into account the presence of a magnetic 
field, which is necessary in order to conserve the polarization of the fuel constituents in a 
magnetic fusion reactor 

For very small colliding energies the analysis of polarization phenomena for the reaction 
d + ^ n + ^if e can be carried out in a general form. In the framework of a formalism, 
based on the polarized structure functions, we will point out the observables which have 
to be measured in order to have a full reconstruction of the spin structure of the threshold 
amplitudes. Data on cross section and tensor analyzing power exist, at threshold [l^ (for 



a review see ||Tj])- Among the two-spin observables, the measurement of a spin correlation 
coefficient, together with the cross section and the tensor analyzing power, allows to realize 
the complete experiment. 

B. Spin structure of the matrix element 

Let us first establish the spin structure of the matrix element. From the P-invariance 
of the strong interaction and the conservation of the total angular momentum, two partial 
transitions, for d + ^He ^ p + '^He (as well as for d + ^H ^ n + ^He) are allowed: 



where Si is the total spin of the He-system and ^/ is the orbital angular momentum of 
the final proton. The spin structure of the threshold matrix element can be parametrized 
in the form: 

^th = 9s^ ■ D + gd{3k ■ D a ■ k - a ■ D), 

where Xi X2 are the two component spinors of the initial ^He and final p, D is the 
3- vector of the deuteron polarization (more exactly, D is the axial vector due to the positive 
parity of the deuteron) , k is the unit vector along the 3-momenta of the proton (in the CMS 
of the considered reaction) and Qg and are the amplitudes of the S— and D— production 
of the final particles, and they are complex functions of the excitation energy. Note that, 
in the general case, the spin structure of the matrix element, for the considered processes, 
contains six different contributions and the corresponding amplitudes are functions of two 
variables. 

The general parametrization of the differential cross section in terms of the polarizations 
of the coUiding particles (in S-state) is given by: 

^{d + me) = (^] [ 1 + AliQabkakb) + A^S ■ P 



dQ \ dfl J Q 



+ A-ik -Pk- S + aJ -PxQ], Qa^ Qabh, (3) 



— * 

where {da/dQ)o is the differential cross section with unpolarized particles, P is the ax- 
ial vector of the target {'^He) polarization, S and Qah are the vector and tensor deuteron 
polarizations. The density matrix of the deuteron can be written as: 

1.. 3 



DaDb — '^i^ab — -i'^abcSc — Qab), Qaa — 0, Qab — Qba- (4) 



After summing over the final proton polarizations one can find the following expressions: 

(5) 



A(^^j = -2ne gsg*a-\gd\^ A2(^^'^ ^ -\gs\^ -Tie g,g2 + 2\gd\^ 
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The coefficients Ai are related by the following linear relation: Ai -\- A2 -\- A'i — —1 for any 
choice of amplitudes Qs and g^- The integration of the differential cross section over the 
A;- directions gives: 

a{d + m) = ao(l + AS-P), A = A2 + \a3^ 

and it is independent from the tensor deuteron polarization. 

The presence of S-wave contribution (the amplitude gs), decreases the value of the integral 
coefficient A whereas, in the fusion resonance region, where the D-wave dominates, the 
maximum value, A = 1/2, is reached. In the complete experiment (which gives \gs\'^, Is'dP 
and Re gsgd)-i the amplitudes \gs\ and \gd\ can be found in a model independent way, with 
the help of the following formulas: 

da\ 



9|^,|' = (5 + 2A-4A) 1 







9|^,r = (2-A + 2A)(^^) , (6) 

-QUegsgl = (1 + 4A + ^2) • 

One can see that the integral coefficient A can be determined from polarized nuclei collisions 
by measuring: 

• the tensor analyzing power Aiva.d + ^He — > p + *^i?e, 

• the spin correlation coefficient Cxx = Cyy = A2 (if the 2;— axis is along /c-direction.) 

Let us study now the polarization properties of the outgoing nucleons. We will show 
that can be predicted only from the tensor analyzing power, Ai- The polarization Pf of 
the produced nucleon depends on the polarization P of the initial ^He (or ^H) as follows: 
Pf = piP + p2k k ■ P, where the real coefficients Pi, i = 1,2, characterize the spin transfer 
coefficients (from the initial ^He or to the final nucleon): K^' = pi +p2Cos^9, — 



— * — * 

p-2 sva.9 cos9, where 9 is the angle between k and P. Averaging over the polarizations of 
the initial deuteron, we can find: 




alio 

3pi = -1 + 2Ai, P2 = -2Ai, 3pi +P2 = -1- 

This analysis holds in the presence of S-state only, in the entrance channel. The validity 
of this assumption can be experimentally verified with the measurement of T-odd one-spin 
polarization observables, as the analyzing powers in d + ^He p + ^He induced by vector 
deuteron polarization or d + ^He —>■ p + '^He. This observable is very sensitive to the presence 
of even a small P-wave contribution, due to its interference with the main amplitude. 



C. Helicity amplitudes 

We calculate here the helicity amplitudes Fx^x^Mi with Ai = A^, A2 = X^hci -^3 = \ (or 
A„), in terms of the partial amplitudes gg and g^- This formalism is very well adapted for 
the analysis of angular distributions of the reaction products, in conditions of fusion reac- 
tors (with polarized fuel) and to the description of polarization phenomena. The direction 
of magnetic field B can be chosen as the most preferable quantization axis (2;— axis). The 
formalism of the helicity amplitudes allows to study the angular dependence of the polar- 
ization observables, relative to B. For example, the polarization properties of the neutron 
m. d -\- ^ n -\- ^He can be easily calculated in terms of these amplitudes. 

The peculiar strong angular dependence of all observables is due to the presence (in 
conditions of fusion polarized reactor) of two independent physical directions, k and B. So 
even for the S-state interaction, a non trivial angular dependence of the reaction products 
appears, i.e. some angular anisotropy, related to the initial polarizations. As all the po- 
larizations of both colliding particles depend on the same magnetic field B, the results for 
these observables depend only on the angle 9, between k and B. The case of the collision 
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of polarized beam with polarized target, where the beam and the target may have different 
directions of polarization is more complicated, but it can also be treated in the framework 
of the helicity formalism. 

The deuteron polarization vector D'^^^ (with a definite helicity A), can be chosen as : 
L>(°) - (0,0, 1) and L>W = l/v^(±l, i, 0). So the following expressions for the six possible 
helicity amplitudes can be found: 

3 

Fo+,+ = - (1 - 3cos2 9)gd, = sin^ Oga, 

Fo+,_ = ^ sin 20gd, = ^ sin 20gd, (7) 

= sin 2egd, = [2^, + (1 - 3 cos^ e)gd] , 

where 9 is the nucleon production angle, relative to the B direction. Other possible helicity 
amplitudes, with reversed helicities of all particles, can be obtained from (7), by parity 
reversion. 

One can see, that for coUinear kinematics, i.e. for 9 = 0°, only two helicity amplitudes 
are nonzero: Fo+,+ — 9s + '^gd, ^-+,- — ~V^{9s~ gd)- So for the corresponding differential 
cross section (for collisions of unpolarized particles) one can find: 

The same formula is also correct in the general case (i.e., for 9 0) for collisions of unpo- 
larized particles: 

( ^\ = 1 V 17^ \^ = ( —\ 

This relation between the differential cross section for coUinear kinematics and the cross 
section at any angle 9 is valid for S-state interaction, only. 

D. Collision of polarized particles 

The angular dependence of the reaction products in d + '^H n + "^He for different 
polarization states of the colliding particles can be derived from (7). 
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• Collisions of longitudinally polarized deuterons (A^ = 0), with polarized or ^He: 
ao+(e) = |Fo+,+ |2 + |Fo+,_|2 = Ig.f + 2Re g,g*(-l + Scos^ 6) + \ga\\l + Scos^ 6). (8) 

• d + ^He collisions with parallel polarizations (relative to B): 

a++(^) = + = ^l^.l^sin^^. (9) 

• d + ^He collisions with antiparallel polarizations: 

= + = 2\g,\' + 2Re g,g*il - 3 cos' 0) + ^(1 + Scos^ 0)\ga\''. (10) 

The sum of all these polarized cross sections is independent from polar angle 0: the unpo- 
larized cross section is isotropic, as expected for 5"— state interaction. 

For the pure fusion resonance (with gs = 0), the angular distribution of the reaction 
products depends specifically on the direction of the polarizations of the colliding particles: 
the sin^ ^-dependence for parallel (++) coUisions, becomes a dependence in (1 -|-3cos^^) 
for (+- ) and (0+) collisions, to be compared with the isotropic behavior of the unpolarized 
collisions. Such definite and strong anisotropy can play a very important role in the design of 
the neutron shield of a reactor and of the blanket, where energetic neutrons (from d + ^H ^ 
n + '^He ) can produce through the reaction n + %i ^ + '^He. Once a d + ^i7-reactor 
is beginning to operate, ^if-fuel can be produced in ^Li-blanket. In principle, this blanket 

— * — * 

can contain polarized ^Li, for a more efficient i?-production in n -|- ^Li-collisions. 

From Figs. 1 and 2, one can see that the angular dependence of the cross sections for 
polarized collisions, is essentially influenced by the presence of the S-wave amplitude and its 
relative phase. 

Let us calculate now the following ratios: 

^ J+lax,xMdcos0 
Cidcos0{da/dn)o, 

from Eqs. (8-10) for (TaiA2(^): 
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0+ 



1, R4 



2(4-3/). 



(11) 



So we can write the following limits: 

0<R++< 3/2 (gs = 0), 1/2 <R+-<2 {g^ = 0). 

In the fusion resonance region, {f=l)^ the (++)-conisions increase the reaction yield (in 
comparison with collisions of unpolarized particles) with a maximum coefficient < 3/2, for 
pure D-wave fusion resonance. Using the notations of |jl| one can obtain the following general 
formula for the differential cross section of d + (or d + 3iJe)-collisions: 



c 

+I2 



7-' 



3 cos 



J - (1 -3cos^^) 
( iRe gsg*d 



1 + 



2i?e gsgg 



+ 



(12) 



Here a = d+t+ + d^t^, b = do, c = d+t+ + d_t+ and d+, do, d_ are the fractions of deuterons 
with polarization respectively parallel, transverse, antiparallel to B , while ^4. and t_ are 
the corresponding fractions for . The relations (i+ + rfo + c?- = 1 and t+ + 1_ = 1 hold. The 
case a = 6 = c= 1/3 corresponds to unpolarized collisions. 

Note that the predicted angular dependence for h and c contributions, Eq. (12), differs 
essentially from the corresponding expression of [|l|. It coincides only for the special case 
/ = 1, (^s = 0. The denominator in Eq. (2) from ref. must be also different. 

From (12) one can find the following expression for the differential cross section of colli- 
sions of polarized deuterons with unpolarized ^H: 



da 

dn 



id + 'H) = 2\gd\ 



1 , ^^^ l-3cos^g ^ 2Reg,g*, 



J "^^ 4 V \9d\ 
i.e. it depends on the tensor deuteron polarization only. We used above the standard 
definition: P^z = d+ — 2do + c?_. Due to the (1 — 3 cos^ 6) dependence, after integration over 
9, the cross section, again, does not depend on the deuteron polarization. 



-^In particular the ratio of amplitudes / = 2|5f(ip/(|(7sp + 2|g'^p) was firstly defined in Q. 
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E. Polarization of neutrons in d + collisions 



Using the helicity amplitudes (7) it is possible to predict also the angular dependence 
of the neutron polarization in d + ^ n + ^He, in the general case of polarized particle 
collisions: 

(n+ - n_)^{d + m) = -{d_t_ - d+t+) sin^ ^(1 - 2 cos^ e)\gd\^+ 



+doit+-t_) 1^^1^-2(1-3 cos^ 9) Re g^gl + (1 - 15 cos^ ^ + 18 cos^ e)\gd\ 



+ 



- d-U)- [4|^/ + 4(1 - 3cos^^)ii:e gsg*d + {l - 15cos^^ + 18cos^^)|^d|^ 

where n± is the fraction of neutrons, polarized parallel and antiparallel to the direction of 
the magnetic field. 

Let us write some limiting cases of this general formula: 

(a) Collisions of polarized deuterons with unpolarized ^i?-nuclei: 

da , T* 



K - ^-)^(^ + = (^+ - [l^-l' + (1 - 3cos' g,gl- 



(2-3cos2e)|^d| 



(13) 



(b) Collisions of unpolarized deuterons with polarized '^if-nuclei: 

|5,|2 + 4(l-3cos'e)i?e^,^^ 



+2(2-3cos2 9)\gd\ 



(14) 



In the case of fusion resonance {gg — 0), these formulas reduce to: 

da ^ ^ 9 
(n+ - n_) — (d + sin^ ^(1-2 cos^ 9) {d+t- - d+t+)+ 

dil 4 



do{t+ - t-) + -{d+t- - d-t+) (1 - 15 cos^ ^ + 18 cos^ 9) 



(15) 



Averaging over the polarizations of d (or ^i7) one can find particular expressions: 
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and 

(n+-n_)^(rfV/J) = --{d_-d+){2-3cos^e). 

The angular dependence of most of these polarization observables is sensitive to the relative 
value of the gs and gd amplitudes, due to the (yf^^'^-interference contributions. Of course, in 
the region of the fusion resonance the gd amplitude is dominant. However the temperature 
conditions, typical for a fusion reactor, correspond to collision energies lower than the energy 
of the fusion resonance. Even a small g.s/gd ratio can change the angular behaviour of the 
polarization observables. In Figs. 3-4 we show, in a 3-dimensional plot, the dependence of 
the neutron polarization on the ratio x = \gs\/gd\ and on the production angle 6 for three 
values of the relative phase 6,6 = 0, tt/2, tt, for d and d + ^if-collisions. 

The exact determination of the parameters x and 6, is crucial for thermonuclear processes. 
This is a reason to perform a complete experiment for this reaction as discussed earlier [§. 
The important point is that even at very low energies, where the spin structure is simplified, 
a complete experiment must include the scattering of polarized beam on polarized target. 
The full reconstruction of the threshold matrix elements requires this type of experiment. 



III. PROCESSES d + d^n +^He AND d + d^ p 

A. Introductory remarks 

The d + d —>■ n +^ He and d + d —>■ p +^ H reactions at low energy have a very wide 
spectrum of fundamental and practical applications, from the discovery of tritium and helium 



isotopes [|I5|, to the important role for primordial nucleosynthesis in the early Universe and 
fusion energy production with polarized and unpolarized fuel |T],|TB|. These processes are of 
large interest in nuclear theory: for example, in a four nucleon system, contrary to three 
nucleon system, broad resonant states can be excited [|l9l. The angular dependence of the 



differential cross sections [|20| , ^ and the polarization observables [11-14] for these charge 



symmetric reactions constitutes a good test of the isotopic invariance for the low energy 
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nuclear interaction. The dd— interaction is also connected to muon catalyzed processes 
[lidd) ^ fi + p or (fidd) —^ii + n +^He p6|, where only the P-state of the cici-system 



is present, at low energy. 

In the general case the spin structure of the matrix element for d + d n +^He [p +^H) 
is quite complicated, with 18 independent spin combinations, and therefore with 18 complex 
scalar amplitudes, which are functions of the excitation energy and the scattering angle. 
However, at very small collision energies, where the S*— state deuteron interaction has to 
dominate, this structure is largely simplified. The identity of the colliding deuterons, which 
are bosons, is an important guide for the partial amplitude analysis in order to determine the 
spin structure of the reaction amplitude. The determination of the polarization observables 
is indispensable, for this purpose. The four possible analyzing powers for (i + (i-collisions. Ay, 
A-zz, and ^2,2, — Ayy were measured at < 100 keV, as well as the angular dependence 



of the differential cross section 20,21,25 



The knowledge of the relative role of different orbital angular momenta (and the corre- 
sponding partial amplitudes) is essential for the solution of different fundamental problems 
concerning these processes, like the possibility to build a thermonuclear "clean" reactor with 
polarized d +^ife-fuel. The main reaction d +^He — > p +^He does not produce radioactive 
nuclei, and the possibility to decrease the cross section oi d + (i-collisions (which produces 
n +^He or p +^H) with parallel polarizations, will decrease the production of neutrons and 
the tritium. Direct experimental data about d + d- low energy collisions are absent, so the 
dependence of the cross section on the polarization states of the colliding particles can be 
calculated only from theoretical predictions or from different multipole analysis. 

The theoretical predictions and the results of multipole analysis seem very controversial 
now, even at very low energy. In the first partial wave analysis [0,^] it was found that 
the S-state c/c/— interaction in the quintet state (i.e. with total spin Si = 2) is smaller in 
comparison with the Si = interaction. This was consistent with the conclusion of ref. [Q, 
that in a polarized reactor it is possible to suppress c? + (i-collisions. Later p3, it was pointed 



out that strong central forces with D— state in ^He can induce a large dd— interaction in the 
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quintet state and resonating-group calculations found that polarized collisions are not 
suppressed. On the other hand, DWBA calculations give a large suppression for the ratio 
of polarized on unpolarized cross section, aj^j^/ao 0.08 in the range Ed = 20 — 150 keV, 



even after inclusion of the ^He D-state. A more recent analysis ^T],^ based on i?— matrix 
approach, concludes that this ratio does not decrease with energy. Note that in principle, it 
can be energy dependent [^]. 

Again, a direct measurement of polarized drf-collisions would greatly help in solving 
these problems and the complete experiment will allow to reconstruct the spin structure 
of the reaction amplitude. Therefore, the considerations based on 5*— wave only, have to 
be considered as the first necessary step which can illustrate the possible strategy of the 
complete experiment for this case. 

B. Partial amplitudes 

We establish here the spin structure of the threshold matrix element for the d + c? — *■ 
n +^He {p +^H) process. For S-state dci— interaction the following partial transitions are 
allowed: 



= - 


J"" 


= 0+ 




= 0, 


^/ 


to 


= 2 - 




= 2+ 




= 0, 




= 2, 


= 2 - 




= 2+ 




= 1, 




= 2, 



where Si is the total spin of the colliding deuterons, ^/ is the orbital angular momentum 
of the final nucleon. Note that the Bose statistics for identical deuterons allows only even 
values of initial spin, that is = and Si = 2 for the S-state. The resulting spin structure 
of the threshold matrix element can be written as: 

M = t{xlcr2Xi) [giDi ■ D2 + 92{'ik ■ k ■ D2 - ■ 

+g3{a -k X Di k ■ D2 + a ■ k x D2 k ■ Di 



(16) 



15 



where xi Xs ^ire the 2-component spinors of the produced nucleon and ^He (or ^H), 
Di and D2 are the 3-vectors of the deuteron polarization, k is the unit vector along the 
3- momenta of the nucleon (in the CMS of the considered reaction). The amplitudes gi and 
g2 describe the production of the singlet n +^ife-state, and the amplitude g^- the triplet 
state. The complete experiment in 5— state (id-interaction implies the measurement of 5 
different observables, to determine 3 moduli and two relative phases of partial amplitudes. 

The validity of the S-state approximation in the near threshold region can be checked 
by measuring any T-odd polarization observable, the simplest of which are the one-spin 
observables as the vector analyzing power in the reaction d + d ^ n+^He [^]. Note that Eq. 
(16) is correct also for the threshold matrix elements of the inverse process: n+^He ^ d + d 
{01 p+^H d + d). 



C. Helicity amplitudes 

In order to establish the angular dependence of the reaction products, for collisions of 
polarized particles, in the presence of magnetic field, let us derive the helicity amplitudes. 
The spin structure of the d + d reactions is more complex in comparison to d +^He. The 
analysis of polarization phenomena is also more complicated. It was mentioned in [|l|, 
that an enhancement factor, equal to 2 can be obtained in a polarized plasma 0, for the 
reaction d + d ^ n +^He, if the deuterons are polarized transversally to the direction of the 
magnetic field, i.e. for (OO)-collisions, in an ordinary thermal ion distribution. Alternatively, 
if colliding beams or beam and target methods are used (inertial fusion), the two ions should 
be polarized in opposite direction, relatively to the field. In case of collisions of deuterons 

with parallel polarizations i.e (++) or ( ), a large suppression of the reaction rate is 

expected. 

It is then interesting to analyze all possible configurations of the polarization of the col- 



^Note that this holds only for the partial wave analysis [27, 2. 
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liding deuterons. We can classify the helicity amplitudes according to the following scheme: 

I) 00 collisions : the polarization is transverse to the magnetic field 2 independent am- 

plitudes; 

II) +- 1- collisions : the polarization parallel to the magnetic field — > 4 independent ampli- 
tudes; 

III) H — collisions : collisions with deuterons with antiparallel polarization, in the same 
direction as the magnetic field — > 4 independent amplitudes; 

IV) 0+ collisions : collisions of one deuteron with polarization transverse to the magnetic 
field with the other deuteron polarized along the magnetic field 4 independent 
amplitudes; 

The corresponding helicity amplitudes ^AiA2,A3A4) (with Ai = A,^^, A2 = A^j, A3 = Xsh^, 
A4 = Atv) are given in terms of partial amplitudes: 

(I) ^00,++ = - sin 2eg3, ^00,+- = ^1 - (1 - 3 cos^ e)g2, 

3 

(II) = sin 29g3, J^++,+- = sin^ e{-g2 + gs), 

3 

^++,-~ = 0, ^++ -+ = sin^ e{--g2 + 53), 

(III) ^+_,++ = ^+_,__ = sin 29gs, (17) 

^+-+- = -^+--+ = -yi - ^(1 - Scos^ e)g2. 
(IV) ^o+,++ = ^(-l + 3cos^%3, J^o+,+- = ^ sin 2^(3^2 + ^3), 
•^0+,— = sin^ 9g3, J^o+-+ = sin 2^(-352 + ^3)- 

— * 

where 9 is the nucleon production angle relative to B direction. 
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D. Angular dependence for collisions of polarized deuterons 



After summing over the polarization states of the produced particles, the cross section 

— * — * „ 

of the process d + d ^ n + He, for definite deuteron polarizations, can be written as: 

aoo{e) = 2 (|J^oo,++|' + 1^00,+- r) = 2\gi - g2{l - Scos^^)^ + Ssin^ ^cos^ ^j^al', 

■9 



A3, -^4 



-sin2^|^2r + 2(l + cos'e)|^3|' 



a+.{9) = Ea3,A4 \^+-,X3xf = 2|^i|2 + 2Re gig*{l - Scos^ 9 + (18) 
-(1 - 3 cos^ ey\g2\^ + 2 sin^ 6 cos^ e\gs\\ 

ao+{e)= ^ |j^o+,A3A4r = 9sin2^cos=^^|^2r + (l-3cos2^ + 4cos^^)|^3r. 

A3,A4 

With the help of these formulas we can estimate the corresponding integral ratios: 



Jl^l d cos 6{da/d cos 6*) ' 
which characterize the relative role of polarized collisions with respect to unpolarized ones: 

_3 15 + 4r ^ _36 r _1215 + r ^ _9 r 

where r = (3|g'2p + 2|g'3p)/|g'ip). It is interesting that all these ratios depend on a single 
contribution of the moduli of the partial amplitudes, the ratio r > 0. The ratios Rx^x^ are 
limited by: 

1.2 < i?oo < 3, < R++ < 3.6, 1.2 < < 12, < Ro+ < 0.9, 

where the upper limits correspond to = S's = 0, (when only the gi amplitude is present), 
and the lower limits correspond to g^i = (for any amplitudes g2 and g^). But the exact values 
of i?AiA2 depend on the relative value of the partial amplitudes, through one parameter, r. 

— * — * 

The general dependence of the differential cross section for d+d-coUisions, can be written 
in terms of partial cross sections ax^^x^ as follows: 

^{d + d) = {dl + d^_)a++{e) + dlaoo{e) + 2d+d_a+_(^) + 2do{d+ + d-)ao+{e), (20) 
dil 



where we used the evident relations between ax^x^: a++{9) — a (9), (7o+(^) = ao-{0), 

cr+_(6') = (T_+(^), due to the P-invariance of the strong interaction, and the standard nota- 
tion: d+, do and d- for different deuteron fractions in polarized plasma. 

Using Eq. (20) one can find some interesting limiting cases. Setting for example, d+ — d- 
(deuterons with tensor polarization only: P^z — l—3do, — 0), one can obtain the following 
dependence of the differential cross section on P^z'- 

^{d+ d) = ao{e) + 2Pzzai{e) + \plza2{e), (21) 

where the coefficients ai{9),i = — 2, are linear combinations of the helicity cross sections 

Qaoie) = 2 [a++{e) + a+_(^)] + aoo{0) + 4a+o(^), 

9ai{e) = a++{e) + a+4e) - aooie) - a+oi0), (22) 

— * — * 

So, measuring the P^;2-dependence of the cross section for d+d collisions, one can determine 
all 3 coefficients ai{6) (at each angle 9). This allows to determine the individual hehcity 
partial cross sections cr;^iA2(^)- 

(Too(^) = ao(^)-4ai(^) + 2a2(^), 

ao+{0) = ao{e) - a,{e) - a2{e), (23) 

a++{e) + a+_{e) = 2ao{9) + 4ai(^) + a2{9). 

In order to disentangle the a++{9) and a^_{9) contributions, an additional polarization 
observable has to be measured, from the collisions of vector polarized deuterons {d± = 

da ^ ^ 

— (d + d) = aoi9) + -fia^M - ^^-m- (24) 
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The linear Pz contribution is forbidden by the P-invariance of the strong interaction. 
Only the measurement of the contribution allows to separate the cross sections a^+{6) 
and cr+_(6'). 

This analysis is equivalent to the discussion of the complete experiment (in terms of 
helicity cross sections crAiA2(^))- 

Finally let us derive the polarization properties of the neutrons in the process d + d ^ 
n +^He. Using eqs. (17) for the helicity amplitudes, one can find for the 6 dependence of 
the neutron polarization (for the different spin configurations of the colliding deuterons): 

{n+ - n_)a++{e) =2sin^^(i^ [si^e^s^a + 2 cos^ ^l^sT] , 

{n+ - n.)ao+{9) = 2doc^+ cos^ [-(1 - 2 cos^ ^)|^3r + 3 sin^ ^i?e ^2^/3], (25) 

(n+ - n_)aoo(6') = (n+ - n_)a+_(6') = 0, 

where n+ and n_ are the fractions of polarized neutrons with spin parallel and antiparallel 
relative to the B direction. 

The production of unpolarized neutrons for 00-collisions of deuterons results from P- 
invariance, and for — h collisions results from the identity of colliding deuterons and from 
the P-invariance. 



E. Complete experiment for d -\- d — > n -\-^ He 

Due to three complex partial amplitudes for the S-wave dd— interaction for the process 
d + d ^ n +^He, the measurement of a large number of observables is necessary, in order 
to perform the complete experiment. This study will be based on the formalism of the 
polarized structure functions, previously used in for the process d — > n +^He. 

Let us consider the collisions of polarized deuterons d + d ^ n +^He. The differential 
cross section can be parametrized in the following general form: 
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^= (^] \l + Ai{k-Qi + k-Q2)+A2Si-S2 + A3k-Sik-S2 
ail \dil J Q L 

+A4k ■Qik-Q2 + A5Q1 ■ Q2 + A6QlabQ2ab ^^^^ 
+A7{k ■ SiX Q2 + k- S2X Qi)j , Qia = Qlahh, Qla = Q2abkb, 

where Si and S2 {Qiab a-nd Q2ab) a-re the vector (tensor) polarizations of the coUiding 
deuterons. The real coefficient Ai describes the tensor analyzing power in d + d n +^He, 
A2 — A7 are the spin correlation coefficients in d + d ^ n+^He. The coefficients Ai —Aq are 
T-even polarization observables and A7 is the T-odd one (due to the specific correlation of 
the vector polarization of one deuteron and the tensor polarization of the other dcutcron). 
Note that these coefficients Ai can not fix the relative phases of the singlet amplitudes gi 
and g2 (from one side) and the triplet amplitude gs (from the other side). The complete 
experiment has to be more complex than the determination of the polarization observables 
Ai- The polarization transfer coefficients from the initial deuteron to the produced fermion 
(n or ^H) have to be measured, too. 

After summing over the polarizations of the produced particles in d + J*— > n +^He, the 
following expressions can be found, for the coefficients Ai, i = 1 — 7, in terms of the partial 
amplitudes gk, k — 1 — 3: 

A2 = -\gi\^ + 2I52I' + Igsl^ - Re gig*2, 

^^(S„='''''^-"'''^ '''' 

^a(^) =-Q\g2? + QRegig*2 + 2\g3\\ 
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where {da/dQ)^ is the differential cross section with unpolarized particles: 

Using these expressions, the following relations can be found between the coefficients Ai- 

(a) linear: between T-even polarization observables, 

9 17 

A2 + A3 + - A = ^1 + A - - A + - A = 

(b) quadratic, relating the T-odd asymmetry A^ with the T-even coefficients A, i = 1 — 6; 

9 

-{l + Al- A]) =Al + (A + A)^ + 6(AA + A A + A A) 

Therefore, the measurements of {da/dVl)^ and 3 coefficients A, i = 1 — 3, allow to find 
the moduli of all S-wave partial amplitudes g^, A; = 1 — 3, and the relative phase of the 
singlet amplitudes gi and g2: 

da \ 



18|c/i|' = (9- 12A-M 



3 J 



dnj^ 



-18|^2r = (9 + 18A + 6A + 10A) fe) , 



2|^3r = (3 + 3A + 2A + 2A ^ 



18i?e g^gl = {-9Ai + 2A) ■ 

So these measurements can be considered as the first step of the complete experiment for 
the process d-\- d ^ n +^He in the near threshold conditions. 

Using these expressions, one can find the following expression for the ratio r: 

1 + a 2, ^ 

r = 3- — , a = -(3A + A). 

1 — Za 9 

Let us discuss, for completeness, the simplest polarization phenomena for the inverse 
reaction, n +^He d + d, when the deuterons are produced in S-state. In this case this 
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process can be described by the same set of partial amplitudes. The cross section for the 

— * 

collisions of polarized particles, n + ^He d + d, can be parametrized in the following way: 

^{n + me) = j [l + B,P, .P2 + B,k.P^k- P^] , (28) 

where Pi and P2 are the polarizations for n and ^He and Bi, i=l,2, are the spin correla- 
tion coefficients for the inverse reaction, which can be expressed as functions of the partial 
amplitudes: 





\9i 


2 


+ 2 


\92 


2 


3 


9i\ 


^ + Q\92\ 


' + M93 


2 



4 


93 


2 


3 


9i 


2 + 6 


92 


2 + 4 


93 


2 



^1 - ,2 , 12 , 12' ^2 - ol„ 12 , «l„ 12 , A\„ 12- (^9) 



These coefficients are not independent, —B1 + B2 — 1, as a result of the pecuhar spin 
structure of the matrix element for the n +^ife— collision. Comparing Eq. (27) and Eq. 
(29), one can find the following relation: 

^S2 = 3(l + A) + 2(A + ^3). 

In order to find the relative phases of the singlet and triplet amplitudes in d+d —>■ n+^He, 
it is necessary to measure the polarization transfer coefficients (from the initial deuteron to 
the final neutron). The most general parametrization of the neutron polarization can be 
written in the following form: 

P = p^S+ p2k{k ■ S) + pak x g, 

where Pi, i=l-3, are the real structure functions, characterizing the corresponding coefficients 
of polarization transfer: 

4/m(c/i - g2)g^ 

Pi = P2 = 0, Ps = 



3|^l|2 + 6|^2P + 4|^3P' 

Both T-even coefficients pi and p2 are identically zero. This results from the specific spin 
structure of the threshold amplitude, due to the identity of the colliding deuterons and the 
S-wave interaction. The T-odd coefficient ps is sensitive to the relative phases of singlet and 
triplet amplitudes. 
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IV. ^He AND COLLISIONS 



A. Three-particle production 

We discuss here the spin structure and the polarization observables for the following 

processes (at low energy for colliding particles): 

^2n+^He, +^ He ^ n + p +*He 

^He+^He^2p+'^He, +^He^ d+^He 

Due to the presence of identical fermions in initial and final states, the spin structure of 
the symmetric processes, ^He +^ He — >■ 2p +'^He and ^H +^H 2n +^/Je, for the S-state 
interaction of colliding nuclei is built on a single transition: Si — Sf — 0, with the following 
parametrization: 

where xi and X2 (0i and (^2) are the 2-component spinors of the produced nucleons (colliding 
nuclei), go is the singlet-singlet amplitude, which dependence on the energies of the produced 
particles has a dynamical character and can be, in principle, complicated. However a model 

— * — * 

independent expression of the cross section as a function of the polarizations Pi and P2 of 

the colliding nuclei can be given as: —(Pi, P2) = ( -r-\ (I — Pi ■ p2)- This dependence is 

dcu \du! J 

correct for any 3-particle phase-space element dcu, and then, also for the total cross section. 
Therefore the polarization of colliding nuclei decreases the reaction rate, in comparison with 
unpolarized collisions. In a fusion reactor, this property will favour the plasma production 
through the main 0? -|-^i7-reaction, which has a larger Q- value and uses 0?— fuel, preventing 
the waste of the very expensive ^H in 'non effective' ^H +^H collision. 

Due to the non-identity of the colliding and produced fermions, the process ^H +^He — >■ 
n + p He is characterized by two independent transitions, a singlet one, (with ampli- 
tude equal to half of the ^H +^H amphtude, due to the isotopic invariance of the nuclear 
interaction) and a triplet amplitude, then the matrix element can be written as: 
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M^H+^He) = ^go{xy2xMM2h) + ^^i(xkaCr2Xi^)(0i(T2a„(/.2). (30) 

— * — * 

Therefore the dependence of the cross section for + ^ife-coUisions on the polarizations 
of the colhding particles is: 

dw \d^Jo \9o\ +^\gi\ 3 

The coefficient A is determined by the ratio R of the cross sections for and +^He 

collisions (with unpolarized particles): 



I^ll' 


^ + Re gog*i 


l^c 


i|' + 3bi| 


2 



This result, of course, is valid only for the S-state interaction in the considered reactions. 

The relative phase of the partial amplitudes go and gi determines the polarization transfer 
from the initial ^He to the final n: 

Pn = pPi, P = 2 

Note that in the process +^He we neglected the production of final particles in the D- 
state, which is allowed in the general case by the conservation of the total angular momentum 
and P-parity: Si — 1 ^ Sf — 1, £1^2 — 2, where £1 is the orbital angular momentum of 
the np— system and £2 is the orbital angular momentum of the '^He, relative to the center of 
mass of the np-system. The spin structure for these transitions can be written as: 

— * 

where k is the unit vector along the 3-momentum of the He (if £2 — 2, £1 — 0) or along 
the proton 3-momentum (if £1 = 2, £2 = 0). The matrix element for the production of two 
P-waves, £1 = £2 = 1 is comparable, in principle, to the D-wave production. 



B. The deuteron production 

Let us finally discuss the process +^He — > d +'^He, (Q-value=14.3 MeV), which is 
also occurring in a 0? + 0? or 0? -|-^ iife-fusion reactor. The spin structure of the threshold 
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amplitude can be established using the generalized Pauli principle (which is valid at the 
level of the isotopic invariance) . Due to the isospin value / = for the entrance channel, the 
sum Si + ii must be odd, and for S- interaction (£j = 0), we have Si = 1, with the following 
two transitions allowed: Si — 1 ^ J"^ = 1+ — > if — and if — 2, where if is the orbital 
angular momentum of the produced deuteron. Therefore, the matrix element can be written 
as: 

M = [hoa ■ D* + /i2(3a ■ k k ■ D* - a ■ D*)] 02, (31) 

where D is the deuteron polarization 3-vector, ho and /i2 are the partial amplitudes, corre- 
sponding to if — and 2. These amplitudes are complex functions of the excitation energy. 
The complete experiment, for this process, must contain the measurement of at least three 
different observables. 

The dependence of the differential cross section on the polarizations of the colliding 
nuclei, have the same structure as eq. (28), with: 

^'[^) =l\ho\' + ^Rehoh* + 2\h2\', ^2 = -Sdh^l' + 2Re hoh*), 

This results from the absence of a singlet state of the colliding nuclei, in the discussed 
process. 

— * 

After k- integration, one can find: 

a{Pi, P2) = (7o(l + SPi ■P2), B^Bi + ^B2 = ^, 

independently on the relative value of the partial amplitudes. Therefore the dependence of 
the total cross section from the polarizations can be predicted exactly, due to the specific 
spin structure of the threshold matrix element. 

The information that can be obtained from the tensor analyzing power in the inverse 
reaction, d-\-^He — -\-^He, does not give any new constraint on the partial amplitudes. 
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h2 


2 


+ 2Re hoh*2 


1 


^0 


2 + 2|/l2 


2 



The cross section is written, in terms of the tensor analyzing power Ad' 
and 

Ad 

i. e. Ad — -zB2- This means that the measurement of the tensor analyzing power Ad allows 
o 

to define both correlation coefficients, for the direct reaction + ^He — > d + He. 

The complete experiment has to include the measurement of polarization transfer coef- 
ficients in "^H + ^He d + He, where the vector deuteron polarization S can be written in 
the following general form: 

S = siP + S2kk-P, 
with the following expressions for the real coefficients Sj, i — 1,2 : 



\ho\ 


^ - Re hoh2 


\k 


P + 2|/l2 


2 



—Re hQh*2 + 


h2? 




ho 


2 + 2|/l2 


2 



*1 — I, n , .^1, in ) ^2 



So, finally, the complete experiment for ^H+^He — > d+'^He has to include the measurements 
of {da/dfl)o, Bi and Si: 

9|/ior = 2(35i + 4si-2) l^^j , 



-9|/i2r = (35i + 4si -5) 



fda\ 
Wo 



-9Re hoh*2 = 2(-35i - Si + 2) 

The results of all other polarization experiments can be expressed in terms of these quantities. 
As an example, let us consider the tensor deuteron polarization from ^H + ^He —>■ d +*He: 



Qab = Q ka{k X P)h + kb{k X P)a , Q 



Im /io^2 



|/io|' + 2|/i2P' 

This T-odd polarization observable can be connected with the T-even coefficients si and Bi 
through the quadratic relation: 

= - 1) - 2(1 - 51 - B^f. 
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In order to analyze the angular dependence of the reaction products for ^H+^He — > d+'^He in 
a magnetic field B, one can extend the formalism of helicity amplitudes, as for the previously 
studied reactions. 

V. CONCLUSIONS 

We have studied nuclear reactions, at low energies, involved in fusion reactors, with 
special emphasis on their strong dependence on the polarization states of the colliding nuclei. 
The results obtained here on the angular dependence and the reaction rate dependence on 
the polarizations, can be used as a guideline in the conception of magnetic fusion reactors. 
The polarization of the produced particles is also important, as it can help the fusion process 
in a working reactor. For example, in a a reactor based on d +^H-iuel, the intensive flux of 
14 MeV neutrons can be used in the Li— blanket, not only for its heating, with consequent 
production of electric power, but also to produce extra ^H- fuel, through the processes: 
n+^Li -^^ H +^He and n+'^Li — > n+^H +^He. Due to the definite polarization properties 
of these reactions, one can increase, in principle, the yield of ^H. 

We showed that the polarization and the angular distribution of the neutrons, produced 
in the process d+^H ^ n +'^He depends strongly on the relative value of the two possible 
partial amplitudes. The presence of a contribution (even relatively small) of the J"^ = 1/2+ 
amplitude is very important for polarization phenomena. 

For the reaction d + d ^ n +^He (with three independent threshold partial amplitudes) 
the situation is more complicated. The d + d-reactions produce energetic neutrons and 
tritium, and should be suppressed in a o? +^He reactor. 

The detailed information about partial amplitudes of different reactions can be obtained, 
in a model independent way, through the realization of the complete experiment. Even 
at low energy, where the spin structure of all matrix elements is highly simplified, the 
complete experiment includes the scattering of a polarized beam on a polarized target. 
These experiments, which are absent up to now, allow the full reconstruction of the spin 
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structure of the threshold amphtudes. 

The main results contained in this paper can be summarized as follows: 

• We give a model independent parametrization of the spin structure of the threshold 

matrix elements for the following reactions: d + d ^ n +^He, d +^ H n +^He, 
^H+^H ^ n + n +^He, ^He +^He ^ p + p +^He, +^He ^ p + n +^He, and 
^H+^He^d+'^He. 

• The hehcity amplitudes for the processes d + d ^ n +^ He and d +^ H — > n +^ He 
are calculated in terms of partial threshold amplitudes, choosing the direction of the 
polarized magnetic field as the quantization axis. 

• The angular distribution of the reaction products for d + cTand d + -^^e-coUisions shows 
a strong dependence on the polarization of the colliding particles, and it can be very 
important to optimize the blanket and the shielding of a reactor. 

• The polarization properties of neutrons, produced in the processes d+^H — >n +^He 
and d + d^n +^He are derived for collisions of polarized particles. 

• The spin structure of the matrix elements and polarization properties are derived 
forsome A=3 induced processes. 
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Figure Caption 



Figure 1. Ratio (To+/(7oo, as a function of x — \gs\/\gd\ and y — cos 9 for the reaction 

d + '^He — > n +*He,ioi different values of the phase 5: (a) 5 = 0, (b) (5 = | and (c) 5 = tt, 
from Eq. (8). 

— * — * . 

Figure 2. Ratio (T+_/(Too as a function of x and y for the reaction d-\- ^He — > n + He, for 
different values of the phase 5: (a) 5 — 0, (h) S — ^ and (c) 6 — n, from Eq. (10). 

Figure 3. Neutron polarization P„ = [— 2 + 3?/^ + a;cos(5(l — 3?/^) + x^] /(2 + X ) , as a 
function of x and y in d +^i?e-collisions, for different values of the phase S: {a) S — 0, (b) 
S — J and (c) S — n. 

Figure 4. Neutron polarization P„ = |^2(2 — 3y'^) + 4xcos5(l — 3y^) + x'^ /(2 + x^), as a 
function of x and y in unpolarized d +^i7e-collisions, for different values of the phase 6: (a) 
(5 = 0, (b) (5 = f and (c) S^n. 
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